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Abstract. Luminous blue variables (LBVs) are characterized by semi-periodic episodes of en- 
hanced mass-loss, or outburst. The cause of these outbursts has thus far been a mystery. One 
explanation is that they are initiated by K-effect pulsations in the atmosphere caused by an 
increase in luminosity at temperatures near the so-called "iron bump" (T ~ 200,000 K), where 
the Fe opacity suddenly increases. Due to a lag in the onset of convection, this luminosity can 
build until it exceeds the Eddington limit locally, seeding pulsations and possibly driving some 
mass from the star. We present some preliminary results from a parameter study focusing on 
the conditions necessary to trigger normal S-Dor type (as opposed to extreme jy-Car type) out- 
bursts. We find that as Y increases or Z decreases, the pulsational amplitude decreases and 
outburst-like behavior, indicated by a large, sudden increase in photospheric velocity, becomes 
likes likely. 
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1. Introduction 

Luminous blue variables (LBV) represent a short-lived phase of massive star evolution 
characterized by episodes of intense mass loss, or outburst. There are two main types of 
outburst behavior. The rare eruptions such as happened with ry-Car and P-Cyg, occur 
every few decades to centuries. They last as long as a few tens of years and blow off ^ IMq 
of material in a single episode. The star also experiences an increase in luminosity during 
this eruption phase. 

In contrast, the more normal S-Dor-type episodes occur every few months to few years. 
They eject ~ 10 _ Mq per event. The LBV does not become more luminous during these 
episodes. For the purposes of this paper, we will call the more extreme 77-Car-like episodes 
eruptions and the less extreme S-Dor-like episodes outbursts. 

Vink & de Koter (2002) determined that the outburst phase results from changes 
in the dominant ionization states of Fe in the wind which result from changes in the 
stellar radius and effective temperature, but what causes the changes in the radius and 
temperature is not clear. In this paper we explore the hypothesis that it is due to seeding 
from K-effect pulsations arising from a large increase in luminosity due to the increase in 
Fe opacity from the so-called "hot iron bump" at temperatures near 200,000 K. We will 
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Table 1. Properties of analyzed modes 



Y 


Z 


Mode 1 


Fractional Energy 
Gain / Period 


29 


0.01 


1H 


i.ii 


29 


0.02 


1H 


2.07 


49 


0.01 


1H 


0.552 


49 


0.02 


1H 


2.22 



also discuss the possiblity that a super-Eddington luminosity in regions of the star near 
these iron bump temperatures could directly drive mass from the star. 



2. Models 

We start with massive main sequence models in the mass range 50M Q ^ M* ^ 80M Q . 
We evolve the stars using the Iben evolution code (Iben 1965). Important additions to 
the code are the Swenson SRIFF analytical EOS, OPAL opacities using the Grevesse 
and Noels 93 mixture (as used in Iglesias & Rogers (1996)), and the empirical massdoss 
prescription of Nieuwcnhuijzen & de Jager (1990). We vary the massdoss rate to allow 
different surface abundances at the end of the evolution. More details of this and the 
other codes used in this analysis are found in Guzik, Watson, & Cox (2005). 

Once we have a model with properties similar to an LBV star, we run the model 
through a pulsation analysis code. This code determines the first three harmonics of the 
K-effect pulsations. The most unstable mode - that is, the one with the largest growth 
rate, as determined by the fractional energy gain per period - is then chosen for further 
analysis. The modes chosen for this analysis are shown in Table [TJ We analyze this 
unstable mode using the ID Lagrangian hydrodynamics code Dynstar (Cox & Ostlie 
1993). This code includes the time-dependent convection treatment formulated by Ostlie 
(1990) and described below. 

The models have log L/L G » 6.1 and T eff = 1.6 x 10 4 K. We vary the metallicity 
Z and the helium abundance Y. Since the pulsations are driven by iron opacity, they 
are Z-dependent, with lower Z resulting in smaller amplitudes. As Y is increased, the 
star effectively becomes older and closer to the stable Wolf-Rayet phase, also resultin in 
smaller pulsation amplitudes. 



3. Time-Dependent Convection 

Our time-dependent convection model is described in detail in Ostlie (1990). Here are 
the main highlights. 

The time-dependent convection model is built upon standard mixing length theory 
(Bohm ??; Cox & Giuli 1968). In the standard theory the convective luminosity is 

g£Q 

where r,T,Cp, p, and g are the local radius, temperature, specific heat, density, and 
gravitational acceleration, respectively. The mixing length I is assumed to be equal to 
the pressure scale height. The parameter Q defined 

< 3 - 2 > 
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Q = 1 for an ideal gas with a constant mean molecular weight. The convective eddy 
velocity in the static case is 




(3.3) 



where Hp is the pressure scale height, 



f= h/l + 4A 2 (A-A od )-l]/A, (3.4) 

A — A 0( j is the superadiabatic gradient, and k is the local opacity. 

The static case, however, does not account for the effect of convective lag due to the 
inertia of the material. This lag is included in our model using a quadratic Lagrange 
interpolation polynomial of v c over the current and previous two time steps. 

For zone i and time step n the time-dependent value of v c is 

„ _ (t'-t n ^)(t'-t^ 2 ) (t> - t n )(t> - t n _ 2 ) ra _ x 

(in — t n ^i)(t n — t n -2) ' \t n —x—t n )[t n -i—t n —2) 

(t' - t n )(t' - t n -i) ^„_ 2 
(in— 2 — t n )(t n -2 — t n -l) C '' % 

where u° ^ is the value determined from equation 13.31 

t' = t„-i + r(t„ - tn-i) (3.8) 

and 

_ _ v ci(tn — t n -l) 1 (o Q\ 

T= Ifac- (3.9) 

The lag factor lf ac = l in our models. 



4. Results 

Figure [1] shows how the photospheric velocity v p tot is affected by changing the metallic- 
ity The figure shows the photospheric velocity as a function of time for a model at solar 
metallicity (defined here as Z = 0.02) and at half-solar, or LMC-like Z . Both models 
exhibit "outburst" -like behavior, in that extreme super-Eddington zones deep in the star 
near 200,000 K (see Figured]) result in a large, sudden increase in the surface velocity. As 
expected the model with the larger metallicity experienced a larger velocity amplitude. 

The Z = 0.01 model with Y increased to 0.49 is shown in Figure [31 In contrast to 
the "outburst" models, this model settles into steady pulsations, and Figure 2] shows 
that this model only becomes super-Eddington by a few percent. Thus the star seems 
able to recover and settle into steady, albeit very large amplitude (v p hot — 100 km s~l) 
pulsations. 



5. Conclusions and Future Work 



We have preliminary results from a study on K-effect pulsations in LBV stars, with 
possible implications for S-Dor-like mass loss. We have found that regions in the star 
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Figure 1. Photospheric velocity as a function of time. Solid line is Z — 0.02, dashed line is Z 

= 0.01. Negative velocity is outward. 




near T <~ 200,000 K experience large super-Eddington luminosities due to the increase of 
radiation from iron bump opacities. We have found that this can cause large amplitude 
pulsations at the surface. At lower Y and higher Z, a large sudden increase in radial 
velocity is seen at the surface. Whether the pulsations or the sudden velocity increases 
evolve into S-Dor-like outbursts is not clear. 

The results presented here are still very preliminary. Much more work is needed to 
confirm our results. Not the least of which is a zoning study. We have used 60-zone 
models to improve the throughput, but it is not clear that this is sufficient. Runs at 120 
zones and higher should be run in order to confirm that the results are robust. 

In addition it is not clear that our "outbursts" result in much, if any mass loss directly. 
As Owocki, et al. showed in these proceedings, the most likely scenario for mass loss in a 
super-Eddington atmosphere is one in which the material separates into clumps with the 
clumps essentially stationary and some small amount of material driven above the escape 
speed between these clumps. As we have only a ID analysis thus far, it is impossible for 
us to determine what the clump size distribution would be. 
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Figure 4. Radiative luminosity at T ~ 210,000 K relative to Eddington for 

Y = 0.49, Z = 0.01. 

This is the beginning of a more formal parameter study on the nature of these pulsa- 
tions and their connections, if any, to S-Dor behavior. In addition to more coverage of 
the parameter space, we eventually plan to extend the hydrodynamic analysis to 2D to 
determine whether the super-Eddington nature of the pulsations leads to enhanced mass 
loss through pores in the stellar atmosphere. 
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Discussion 

Langer: How many pulsation cycles does it take to lose the pulsating layer in a station- 
ary wind, assuming standard radiation-driven wind mass-loss rates? 

Onifer: The pulsation cycle for the model I showed was about 1 month, and the Fe 
bump region is <~ 1O~ 3 M0 below the surface, so it would take about 100 years or ~ 1000 
pulsation cycles to strip away the layer. 

Humphreys: Thanks for making the distinction between the normal or classical LBVs 
like S Dor & AG Car and the r\ Car-like giant eruptions. r\ Car & its rare relatives actu- 
ally increase their luminosity during their giant eruptions. The classical LBVs maintain 
constant luminosity during their optically thick wind stage. Owocki was talking about 
?7-Car-like eruptions. It may only be a difference of degree, and the mechanism may be 
the same, but what we observe is very different. It is important to make the distinction. 

Onifer: Thank you, this is a very good point. 

Hirschi: If pulsations are driven by iron opacity, pulsations would be strongly metallicity 
dependent. Are there other sources of opacity at very low metallicities to drive pulsations. 

Onifer: I'm not aware of any sources of opacity at similar temperatures that could take 
over for iron at low metallicity. 



